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The compatibility of the pure power law infrared solution of QCD Dyson-Schwinger equations 
(DSE) and lattice data for the gluon and ghost propagators in Landau gauge is discussed. For the 
gluon propagator, the lattice data is compatible with the DSE infrared solution with an exponent 
K ^ 0.53, measured using a technique that suppresses finite volume effects and allows to model 
these corrections to the lattice data. For the ghost propagator, the lattice data does not seem to 
follow the infrared DSE power law solution. 
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The problem of quark and gluon confinement in QCD is still not understood In the Lan- 
dau gauge, two proposed confinement mechanisms, the Kugo-Ojima and the Gribov-Zwanziger 
scenarios, relate gluon confinement to the infrared properties of gluon and ghost propagators. 

The investigation of the infrared behaviour of gluon and ghost propagators requires a non- 
perturbative formulation of QCD such as the Dyson-Schwinger equations (DSE) or the lattice for- 
mulation of QCD. The nature of the lattice QCD approach allows for a study of the propagators 
including all non-perturbative physics. However, the lattice is limited to a finite volume, therefore 
the access to the deep infrared regime of QCD requires a proper analysis of finite volume effects. 
On the other hand, the DSE (for a recent review see [Q]) were solved analitically in the infrared, 
predicting a pure power law behaviour for the gluon and ghost dressing functions. 



2\2k 



Zghost{q )=q G{q)r^{q) 



2\-K 



(1) 



The solution requires a truncation of an infinite tower of equations and a parametrization of a 
number of vertices. Assuming ghost dominance and a bare ghost-gluon vertex, DSE estimated an 
exponent k = 0.595 This value for K, being above 0.5, implies a vanishing gluon propagator 
and a divergent ghost dressing function, in agreement with the confinement criteria referred above. 
Other studies |^ ||, |6|, ^ |8]] gave further support to this picture. 

On the lattice, in order to test reliably for the power law solutions (|l|) one needs very large 
volumes. A possible cheaper solution could be the use of large asymmetric lattices, i.e. lattices as 
X Lf, with a large Lf — see [Q] and references there in. Although there are non-negligible finite 
volume effects due to the small spatial lattice extension, the large temporal lattice extension allow 
the access to momenta below 100-200 MeV, and a direct test on the validity of the solution (|l]). 

So far, we have computed K" by a direct fit of (|l|) to the asymmetric lattice data. The results 
show that K increases with the spatial lattice volume [^. In this proceeding we report on a technique 



niOD that seems to provide a volume independent value for the gluon exponent K. 

The lattice setup used in this work is summarized in table [I[ The configurations were gauge 
fixed to the Landau gauge using a Fourier accelerated Steepest Descent method, starting from 



Lattice 


Update 


Therm. 


Sep. 


# Conf. 


8^ X 256 


70VR-F4HB 


1500 


1000 


80 


10^ X 256 


70VR-F4HB 


1500 


1000 


80 


12^ X 256 


70VR-F4HB 


1500 


1000 


80 


14^ X 256 


70VR-F4HB 


3000 


1000 


128 


16^ X 256 


70VR-F4HB 


3000 


1500 


155 


18^ X 256 


70VR-F4HB 


2000 


1000 


150 


16^ X 128 


70VR-F2HB 


3000 


3000 


164 



Table 1: Lattice setup. All sets of configurations were generated using a combined Monte Carlo sweep of 
overrelaxation (OVR) and heat bath (HB) updates. The number of thermalization (Therm.) and separation 
(Sep.) sweeps refers to combined sweeps. 
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the identity gauge transformation. Then, the gluon propagator was computed using the same 
definitions as in The ghost propagator was computed using a plane-wave source |11]. Statisti- 
cal errors on the propagators were computed using the jackknife method. Otherwise, the statistical 
errors were computed using the bootstrap method with a 68% confidence level. 

On the lattice, the finite volume and discretization effects can be suppressed by defining ratios 
between similar quantities. For example, consider some quantity A{x). Suppose that the lattice 
effects are given by 1 + with 5 ^ 1 . If x' ~ x, one can write 



A(x)(l + 5) A{x) 



A{x) 



i.e. the error on the ratio is of second order in 8. 

Let's take the gluon propagator. If, in the continuum limit, we have Z{q^) 
finite lattice we have, in general. 



(2) 



(0{q^ f-*^, on a 



(3) 



where ^{q) can be viewed as a multiplicative correction to the continuum function. Then, taking 
ratios of the gluon dressing function between consecutive temporal momenta. 



q[n] =qA[n 

and taking logarithms, one gets 

' Ziatt{q^[n+\]) 
Ziatt{q^[n]) 



2 . /Tln\ 
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191 MeV 


238 MeV 


286 MeV 


333 MeV 


381 MeV 


8 


K 


0.526(27) 


0.531(19) 


0.531(13) 


0.522(16) 


0.527(12) 




c 


-0.179(54) 


-0.194(34) 


-0.193(19) 


-0.171(28) 


-0.184(18) 






0.12 


0.11 


0.08 


0.48 


0.54 


10 


K 


0.511(35) 


0.531(25) 


0.525(21) 


0.523(17) 


0.527(16) 




X^/d.o.f. 


0.69 


0.98 


0.74 


0.56 


0.50 


12 


K 


0.509(31) 


0.517(21) 


0.508(18) 


0.521(18) 


0.530(14) 




c 


-0.094(56) 


-0.112(35) 


-0.094(25) 


-0.119(27) 


-0.138(18) 




xVd.o.f. 


0.11 


0.16 


0.33 


0.84 


1.03 


14 


K 


0.536(24) 


0.540(19) 


0.548(16) 


0.545(12) 


0.542(11) 




xVd.o.f. 


0.33 


0.20 


0.39 


0.34 


0.34 


16 


K 


0.539(22) 


0.528(17) 


0.534(12) 


0.536(12) 


0.539(11) 




c 


-0.125(43) 


-0.102(30) 


-0.112(19) 


-0.118(14) 


-0.123(12) 




x'-fd.o.f. 


1.77 


1.24 


0.96 


0.78 


0.68 


18 


K 


0.529(20) 


0.516(16) 


0.523(14) 


0.536(11) 


0.5398(95) 




c 


-0.099(36) 


-0.068(25) 


-0.085(19) 


-0.111(14) 


-0.119(13) 




X'-Id.o.f. 


0.39 


0.77 


0.85 


1.79 


1.58 



Table 2: Linear fits of /?|"""[n] as a function of Rq[n]. 



Defining 



R 



,gluon 



n] = In 



Z,,niq^[n+\]) 
Ziatt{q^[n]) 



Rq[n]=\n 



q^[n+\] 



equation becomes 



lKRq[n\ + C(^). 



(6) 



(7) 



Figure |l| shows these functions for the asymmetric lattices with = 256. For each lattice 
size, the data points define straight lines in the infrared region, i.e. it seems that C(g) is a constant, 
Cicj) = C. This hipothesis can be tested by fitting R^^"""[n] to a linear function of Rq[n]. The 
corresponding values for K and C are reported in table ^ Note that the measured K values are 
stable against variation of the fitting range and spatial lattice size. Furthermore, within one standard 
deviation, k > 0.5 as in the solution of the DSE. 

The results of the linear fit suggest a parametrization of the finite volume effects. From the 
definition of A, it follows that A(g[?i + 1]) = A{q[n\)e'-^ and 



dA{q) Aiq[n+l])-A{q[n]) 



dq 



q[n+l\ -q[n] 



'Kq)- 



In 
oLt 



A{q)A 



(8) 



where A is a constant. The integration of this equation gives 



(9) 
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Qnmx '• 


191 MeV 


238 MeV 


286 MeV 


333 MeV 


381 MeV 


8 


K 


0.526(26) 


0.533(19) 


0.534(11) 


0.523(10) 


0.524(9) 




A{GeV'^) 


-3.75± 1.1 


-4.06(68) 


-4.11(34) 


-3.69(28) 


-3.73(23) 




X'id.o.f. 


0.09 


0.12 


0.08 


0.62 


0.51 


10 


K 


0.511(27) 


0.536(22) 


0.534(17) 


0.531(14) 


0.534(13) 




X^d.o.f. 


0.53 


1.08 


0.73 


0.58 


0.49 


12 


K 


0.508(31) 


0.515(22) 


0.507(15) 


0.520(12) 


0.537(9) 




A(GeV-^) 


-1.9± 1.2 


-2.25(78) 


-1.92(46) 


-2.40(36) 


-2.96(23) 




X^d.o.f. 


0.07 


0.12 


0.24 


0.84 


1.94 


14 


K 


0.538(23) 


0.542(18) 


0.552(14) 


0.551(11) 


0.546(9) 




X^/d.o.f. 


0.24 


0.17 


0.47 


0.36 


0.45 


16 


K 


0.541(22) 


0.532(16) 


0.535(10) 


0.539(9) 


0.543(8) 




A{GeV-^) 


-2.67(84) 


-2.29(54) 


-2.39(31) 


-2.53(24) 


-2.66(18) 




X^/d.o.f. 


1.15 


0.78 


0.55 


0.50 


0.54 


18 


K 


0.529(20) 


0.516(15) 


0.523(12) 


0.539(9) 


0.550(8) 




A{GeV-^) 


-2.05(79) 


-1.50(51) 


-1.75(33) 


-2.31(24) 


-2.66(20) 




xVd.o.f. 


0.28 


0.59 


0.54 


2.14 


2.71 



Table 3: Fits of the gluon dressing function to the pure power law with exponential correction. 



i.e. the lattice dressing function is given by 

ZLatiq^) = C0{q^f%^''- (10) 
note that now (O has absorved the constant Aq. The finite volume effects are summarized by the 



constant A. The results of fitting the lattice data to ( |10[ ) are reported in table ^. The K values in 
tables ^ and |3| are, as expected, essentially the same. Furthermore, the constants A and C are, in 
general, decreasing functions of the (spatial) lattice volume, and should go to zero in the infinite 
volume limit. 

The ratio method provides similar results when used with the 16^ x 128 lattice data. As shown 
in [^, the gluon data from 16^ x 128 and 16^ x 256 lattices are compatible within errors and one 
expects similar values for the constant A for these two lattices. Indeed, the fits to (10), see table ^, 
give essentially the same A. Given the relation between A and C, 

— \ aL, 

aLf 

we expect, as observed, that Ci28 — 2 x €256- Furthermore, from table ^ one can conclude that the 
effects of Gribov copies are not resolved by the statistical precision of our simulation. 

Note that the estimated gluon infrared exponent K ~ 0.53 implies a vanishing gluon propagator 
at zero momentum. Figure ^ shows the bare D{0) as a function of x = I /V . In what concerns the 
infinite volume D{0), a linear (Doo(O) + bx) or quadratic (Doo(O) + bx + cx^) extrapolations give a 
non-zero Doo(O), but a power law (ax*) extrapolation, which implies Doo(O) = 0, is also possible. 
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100 




2e-06 



4e-06 , 
X = l/V 



6e-06 



8e-06 



Figure 2: Values of D{Q) as a function of 1 /V . 



and gives b = 0.10. Curiously, this value is very close to the figure reported in a recent investigation 
of the DSE on a torus b ~ 0.095. 



In what concerns the ghost propagator, see [10] for details, we were not able to see a linear 
behaviour of as a function of Rq[n\. Assuming that the ghost dressing function really 

follows a pure power law in the infrared, the only result one could extract is a kind of lower bound 
for the ghost infrared exponent, K ~ 0.29. 

Finally, a few words about the running coupling defined from the gluon and ghost propagators. 



ghost 



(12) 



For the lattice data strong coupling constant see figure ^ If the DSE predict a finite non-zero value 
for 0:5(0) = 2.972, the lattice data seems to go to zero. This behaviour is in agreement with the 



solution of the Dyson-Schwinger equations on a torus []12|]. According to this study, one should go 
to even larger lattices to become closer to the continuum. This is a very ambitious challenge for 
the next years. 







Ratios 






Modelling 




16^^ xL, 


K 


C 


xVdof 


K 


A{GeV-^) 




L, = 256 


0.539(11) 


-0.123(12) 


0.68 


0.543(8) 


-2.66(18) 


0.54 


L, = 128 [ID] 


0.541(19) 


-0.239(38) 


0.01 


0.542(20) 


-2.56(39) 


0.01 


L, = 128 [CEASD] 


0.539(19) 


-0.234(36) 


0.15 


0.539(18) 


-2.47(36) 


0.10 



Table 4: Results obtained for the lattices with L.s = 16 < 381MeV). For the lattice 16^ x 128, two gauge 
fixing methods were considered. ID stands for a gauge fixing started from the identity gauge transformation, 
and CEASD stands for the gauge fixing method devised in [[l3|], aiming to find the global maximum of Fu [g\. 
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q{GeV) 



Figure 3: Comparison of results for (Xs{q^) using both lattice QCD and DSE. 
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